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Abstract: We study the Higgs sector of the Next to Minimal Supersymmetric Standard 
Model (NMSSM) in hght of the discovery of the SM-hke Higgs boson at the LHC. We 
perform a broad scan over the NMSSM parameter space and identify the regions that 
are consistent with current Higgs search results at colliders. In contrast to the commonly 
studied "decoupling" scenario in the literature where the Minimal Supersymmetric Standard 
Model CP-odd Higgs boson mass is large rriA ^ rnz^ we pay particular attention to the 
light Higgs states in the case when rriA ^ 2m^. The Higgs bosons in the NMSSM, namely 
three CP-even states, two CP-odd states, and two charged Higgs states, could all be rather 
light, near or below the electroweak scale, although the singlet-like states can be heavier. 
The SM-like Higgs boson could be either the lightest CP-even scalar or the second lightest 
CP-even scalar, but is unlikely to be the heaviest scalar. These NMSSM parameter regions 
have unique properties and offer rich phenomenology. The decay branching fractions for the 
SM-like Higgs boson may be modified appreciably. The correlations of ^^/VV and VV/bb 
can be substantially altered. The new Higgs bosons may be readily produced at the LHC 
and may decay to non-standard distinctive final states, most notably a pair of Higgs bosons 
when kinematically accessible. We evaluate the production and decay of the Higgs bosons 
and comment on further searches at the LHC to probe the Higgs sector of the NMSSM. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) remains one of the most appeal- 
ing models that leads to a more complete theory beyond the Standard Model (SM). Its 
simplest extension, the Next to Minimal Supersymmetric Standard Model (NMSSM) [1, 2] 
introduces additional appealing features. Among the most notable is that it provides an 
attractive solution to the SUSY fi problem [3]. Furthermore, it is widely believed that 
the discovery of a Standard Model-like Higgs boson [4, 5] strongly supports the idea of 
weak-scale supersymmetry (SUSY) based on the "naturalness" argument. However, in the 
context of the MSSM, a Higgs mass of rrih ^ 126 GeV still requires a significant degree of 
fine tuning [6-9]. In contrast, the NMSSM largely alleviates the tuning required to achieve 
this rather high mass value [10, 11]. 

The Higgs sector of the MSSM consists of two SU(2)i, doublets. After electroweak 
symmetry breaking (EWSB), there are five physical states left in the spectrum, two CP- 
even states and with m/^o < m^o, one CP-odd state A^, and two charged scalar 
states H^. At tree-level, it is customary to use the mass uia and the ratio of the vacuum 
expectation values tan/3 = Vu/v^ as the free parameters to determine the other masses. 
These masses receive large radiative corrections from the top-stop sector due to the large 
top Yukawa coupling. If we categorize these Higgs bosons according to their couplings to 
the electroweak gauge bosons, there are two distinct regions in the MSSM [12]: 

(i) The "decoupling region": For a relatively heavy {rriA ^ 300 GeV), the lighter 
CP-even state is the SM-hke Higgs and the others A^, and are heavy and 
nearly degenerate. 

(ii) The "non-decoupling region": For ttia ^ rnz, the heavier CP-even Higgs is the 
SM-like Higgs, while and A^ are light and nearly degenerate. The mass of the 
charged Higgs is typically around 140 GeV. 

The decoupling scenario comfortably accommodates the current searches due to the effective 
absence of the non-SM-like Higgs states. In fact, it would be very difficult to observe any 
of the heavy MSSM Higgs bosons at the LHC if tua ^ 400 GeV for a modest value of 
tan/3 < 10 [13-16]. The non-decoupling scenario, on the other hand, would lead to a 
rich LHC phenomenology due to the existence of multiple light Higgs bosons. Although 
this latter scenario would be more tightly constrained by current experiments, it would 
correspondingly have greater predictive power for its phenomenology. 

In the NMSSM, one complex SU(2)i,xU(l)y singlet scalar field is added to the Higgs 
sector. As a result, after the scalar fields acquire vacuum expectation values, one new 
CP-even and one new CP-odd state are added to the MSSM spectrum, resulting in three 
CP-even mass eigenstates (denoted by i^i, 7^2, ^^3), two CP-odd mass eigenstates (Ai, A2), 
plus a pair of charged states (H^). 

The masses of the CP-even scalars can be better understood by considering what 
happens when the singlet is added to the MSSM spectrum. The masses of the MSSM 
Higgs bosons can be in one of two scenarios: the SM-like Higgs of the MSSM can either be 
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Figure 1. Illustration of the effect of adding the singlet to the MSSM CP-even Higgs boson 
spectrum before mixing. 



the lighter eigenstate or the heavier eigenstate, as illustrated in the top row of Fig. 1. After 
adding the singlet scalar, the two panels of the MSSM give rise to six possible scenarios in 
the NMSSM, as illustrated in the lower row of Fig. 1. In reality, the mass eigenstates are 
admixtures of the gauge interaction eigenstates, and thus cannot be labelled as simply as 
in Fig. 1. Nevertheless, these graphs give us an intuitive picture of the result of adding the 
singlet field of the NMSSM. 

Recently, many analyses of the NMSSM have been performed in light of the recent 
Higgs searches at the LHC, focusing on the large ttia region. References [17-19] showed 
the compatibility of the NMSSM with an enhanced 77 rate, while Reference [20] studied 
the stringent flavor and muon g — 2 constraints on the model. Moreover, the NMSSM may 
include many interesting features, that include grand unification of gauge couplings [21], 
naturalness for the Higgs mass [10, 11, 22-24], neutralino Dark Matter [25-27], and possible 
accommodation of multiple nearly degenerate Higgs bosons [28-30]. 

In this paper, by contrast, we consider the NMSSM in the Iow-tua region, the prototype 



-3- 



of which is the non-decouphng scenario of the MSSM. We perform a broad scan over the 
NMSSM parameter space and identify the low-m^ regions that are consistent with current 
Higgs search results at the cohiders, including the discovery of a SM-like Higgs boson. We 
find that the Higgs bosons of the NMSSM, three CP-even states, two CP-odd states, and 
two charged Higgs states, could all be rather light, near or below the electroweak scale 
in our low-m^ scenario, although the singlet-like states can also be heavier. The SM-like 
Higgs boson could be either the lightest scalar or the second lightest scalar, as illustrated in 
panels 1 — 4 of the bottom row of Fig. 1. However, it is extremely difficult to uncover any 
regions corresponding with the scenarios of the last two panels of Fig. 1 where the SM-like 
Higgs boson is the heaviest CP-even state after imposing all the existing collider search 
constraints. 

These low-m^ parameter regions of the NMSSM have unique properties and offer rich 
phenomenology, providing complementary scenarios to the existing literature for the large- 
rriA case, as mentioned above. The production cross section and decay branching fractions 
for the SM-like Higgs boson may be modified appreciably and new Higgs bosons may be 
readily produced at the LHC. We evaluate the production and decay of the Higgs bosons 
in this model and propose further searches at the LHC to probe the Higgs sector of the 
NMSSM. 

The rest of this paper is organized as follows. In Sec. 2, we present a short, self- 
contained introduction to the Higgs sector of the NMSSM. In Sec. 3, we discuss our param- 
eter scanning scheme and the current constraints applied. We then discuss the resulting 
constraints and correlations for the NMSSM parameter space in Sec. 4 for the case that the 
SM-like Higgs is the lightest CP-even scalar (panels 1-2, bottom row of Fig. 1) and in Sec. 5 
when the SM-like Higgs is the second lightest CP-even scalar (panels 3-4, bottom row of 
Fig. 1). In Sec. 6, we consider the basic LHC phenomenology for our results. Finally, we 
summarize and conclude in Sec. 7. 

2 NMSSM Higgs Sector and the Low-m^ Region 

In the NMSSM [1, 2, 31, 32], a new gauge singlet chiral superfield S is added to the MSSM 
Higgs sector resulting in a superpotential of the form 

W^NMSSM = Yuu'^HuQ + YddFHdQ + Yee'HdL + XSHuHd + (2.1) 

with an explicit Z3 symmetry. Additionally, the soft-SUSY breaking Higgs sector of the 
NMSSM is: 

VH,Soft = ml^HlHu + ml^H\Hd + mI\S\^ + (^XAx(HfeHd)S + ^kA^S^ + c.c)j . (2.2) 

After the singlet obtains a vacuum expectation value (VEV) (5) = Vs/\/2^ an effective /i 
term is generated: ji — \vs/\/2^ which solves the so-called //-problem of the MSSM. An 
effective 6-term 6eff — m(^a + f^) is also generated at tree level. 

In this work, we assume a CP-conserving Higgs potential with all the coefficients being 
real. We further take A and Kj to be positive, unless otherwise stated. For the VEVs, we 
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use the convention {H^) = Vu/V^, {Hf) = Vd/V2, with vl^v^^ v'^ ^ (246 GeV)^ and 
tan/3 = Vu/vd- After electroweak symmetry breaking, we are then left with three CP-even 
Higgs states i^i, 7^2, ^^3, two CP-odd Higgs states Ai, ^2, and a pair of charged Higgs 
states 

2.1 Masses 

2.1.1 CP-odd Higgs Bosons 

For the CP-odd Higgs bosons, we define the mixing states 

A^^V2 {lm{H^) sin/3 + lm{H^) cos/3) , As = V2 Im(5). (2.3) 

The relevant parameters of our interest are the diagonal elements of the mass matrix in the 
basis of {Ay^ Ag) as 

2jl f Hi \ 26eff 



sm2/3 V A / sm2/3 



= (^m^sin2/3 + 6^/x^jsin2/3-3^M/^. (2.5) 

The full mass matrix expression can be found in Ref. [31]. In the limit of zero mixing between 
Ay and A^, tua is the mass of the CP-odd Higgs A^^ as in the case of the MSSM. However, 
in the NMSSM, the mass eigenstates are typically a mixture of Ay and A^, resulting in a 
more complicated mass spectrum and parameter dependence. Although ttia is not a mass 
eigenvalue in the NMSSM, it takes the same form in terms of 6eff as in the MSSM [see 
Eq. (2.4)]. We also note that m^^ has the contribution — 3^/xA^. As a result, to obtain 
positive mass squared eigenvalues, the combination jiA,^ can not be too large and positive, 
in particular, for the small ttia region that we consider in this paper. We denote the mass 
eigenstates as Ai and ^2, where rriAi < ^^2- 

2.1.2 Charged Higgs Bosons 

The charged Higgs bosons in the NMSSM have the same definition as in the MSSM, 
but a new contribution to their mass 

= sin /3 + cos /3, m^± = + - \ i^^f- (2-6) 

The extra A-dependent term leads to a reduction of the charged Higgs mass compared to 
its MSSM value. Requiring > gives an upper bound for A as a function of ttia 

V2 



V 



\<—Jm\ + m^^. (2.7) 



The LEP search limit m^± > 80 GeV [33, 34], as well as the bounds from the Tevatron and 
LHC charged Higgs boson searches, strengthen this upper limit even further, depending on 
the value of tan /3. 
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2.1.3 CP-even Higgs Bosons 

The CP-even Higgs sector is much more comphcated compared to that of the MSSM. It is 
advantageous to define the basis as: 

cos/3 sin/3 \ / (Re(ifO) - Vd) \ 
-sm/3 cos/3 ] I (Re(if°) - Vu) j 




S = V2 (Re(5) 



{2.i 



The benefit of using this basis is that the couphngs of to the gauge sector and the 
fermion sector are exactly the same as that of the SM Higgs. On the other hand, Hy does 
not couple to pairs of gauge bosons at all, and its coupling to the up-type quarks (down- 
type quarks and charged leptons) is proportional to ^7^^ (tan/3) with respect to the SM 
values. The singlet, 5, does not couple to either the gauge bosons or the fermions. While 
the mass eigenstates i?i,2,3 (with rriHi < ^/fs — ^Hs) are typically mixtures of /i^, Hy 
and 5, by knowing the fraction of hy^ Hy^ and S in the mass eigenstates, we have a better 
understanding of their interactions with the gauge bosons and fermions. 

The diagonal entries of the mass matrix for the CP-even Higgs bosons in the basis of 
{hy^ Hy^ S) are given by [31] 



mt 


= m| + 




= m\- 






m| 





m| 



{Xvf 
{Xvf-ml 



sin^ 2/3 , 
sin^ 2/3 , 



m\ sin 2/3 



2-ir 



sin 2/3 + —ji 
A 



(2.9) 
(2.10) 
(2.11) 



Note that the combination jIJ^A,^, that appear in m|, also appeared in [see Eq. (2.4)]. 
While jfJ^A,^ could not be too large and positive in order for to be positive, we see 
that it also can not be too large and negative in order for m| to be positive. This term also 
introduces certain correlation between /i and A^, as discussed in Sees. 4 and 5. For large 
m^, we see that m|^^ grows with m^, while rrih^ remains around the electroweak scale. 
The singlet, on the other hand, is determined by a combination of /x, A^, and m^, as well 
as the dimensionless quantities A, and tan/3. 

The first and foremost effect of the introduction of the singlet and its couplings to 
the MSSM Higgs sector is the extra A-term in Eqs. (2.9) and (2.10), which lifts up the 
mass of the SM-like Higgs in particular, for small tan/3, while reducing Trij^^. In the 
MSSM, for the SM-like Higgs to have a mass of approximately 126 GeV typically requires 
the tree-level Higgs mass-squared m^cos^2/3 to be maximized, which prefers large tan/3. 
In the NMSSM, by contrast, the contribution from ^(Ai;)^ sin^ 2/3 results in small values 
of tan/3 being favored, especially for large A. Consequently, the contribution to the Higgs 
mass from stop sector loop corrections can be relaxed. The left-right mixing in the stop 
sector is no longer required to be near maximal (l^^l ^ ^y6MssQMssu in the MSSM). 

The mixture of the singlet with the MSSM Higgs sector, in particular with hy^ could 
further affect the SM-like Higgs mass. If we consider only the hy-S mixing for simplicity, 
when m^^ > m|, the mass eigenvalues for the SM-like Higgs is pushed up after the diago- 
nalization of the 2x2 mass matrix. This is the so-called "push-up" scenario described in 
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the literature. On the other hand, when m^^ < m|, the mass eigenvalue for the SM-like 
Higgs is pushed down due to the mixing, and is, thus, called the "push-down" scenario. 
Such effects have been discussed extensively in the literature [11, 19], considered almost 
exclusively in the limit of rriA ^ rnz^ which decouples the effect of the MSSM non-SM like 
Higgs Hy^ while focusing only on the mixture of hy and S. The low- lying spectrum in such 
cases includes two CP-even Higgs bosons. Hi and 7^2, as a mixture of hy and 5, with either 
Hi or H2 being the 126 GeV SM-like Higgs, corresponding to the push-down or push-up 
scenario, respectively. In this large-m^ scenario, only one CP-odd Higgs Ag might be light, 
while Ay and H^ are heavy and decouple. Both the push-up and push-down scenarios, 
however, suffer from a certain degree of fine-tuning for the NMSSM parameters if the stop 
masses are relatively light and the left-right mixing in the stop sector is not large [11]. 

2.2 Low-m^ Region 

In this paper, we consider the region of the NMSSM with relatively small tua {ttia ^ 2m^). 
In this region, all the MSSM-type Higgs bosons are relatively light, with m^^ and m|^^ 
relatively close to each other. With an appropriate choice of other NMSSM parameters, 
m| and can be light as well. This could lead to potentially large mixing effects in the 
Higgs mass eigenstates, resulting in possible deviations of the SM-like Higgs couplings to 
the gauge boson and fermion sectors. 

The low-m^ region of the MSSM (the so-called "non-decouphng" region) has been 
studied in Refs. [16, 35-40]. It was pointed out that for rriA ^ rnz^ the heavy CP-even 
Higgs, H^ in the usual MSSM notation, is the SM-hke Higgs. On the other hand, for the 
light CP-even Higgs to be SM-like, rriA is typically large: rriA ^ 300 GeV, in the so-called 
"decoupling" region of the MSSM. However, these observations do not necessarily hold in 
the NMSSM, due to the singlet-induced A-term contribution to and TTT/j^ , as well as 
the singlet mixing effects in the mass matrix. 

If we ignore the singlet mixture with hy and Hy for the moment, and study the conse- 
quence of the extra A-term in the 2x2 {hy, Hy) system in the NMSSM, then, to have the 
heavy CP-even MSSM Higgs be SM-like, m^^ > ra^^, requires 

m\ < m|cos4/3 + (Ai;)^sin^2/3. (2.12) 

Fig. 2 shows the lines in the A versus rriA plane when = ^If^' various values of 
tan/3. For regions above the lines, > TTij^^^ heavy CP-even MSSM Higgs is 

SM-like (up to mixing and loop corrections). For regions below the lines, m^^ < Tnj^^, and 
the light CP-even MSSM Higgs is SM-like. All the lines cross the point rriA = rnz and 
A = \/2mz/v ^ 0.5. For small tan/3 ^ 1, large A (above the tan/3 = 1 line) is preferred 
to realize rrih^ > rriH^, while small A gives rise to ruh^ < tuh^. For larger values of tan/3, 
the curve tilts more and more vertically. For tan/3 > 10, the A dependence becomes rather 
weak and the separation of the two regions is governed by the value of mA- tua ^ rnz 
for ruh^ > tuh^ and ttia ^ foi" ^/i^ < '^i^^? which is similar to the usual MSSM 
case. In Fig. 2, we also include the m^± contours as dashed lines, with the shaded area 
indicating the region ruled out by < 0. Taking into account the LEP bound of 
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Figure 2. Lines for m^^ = m'jj^ for different values of tan /3 in the A versus tua plane, m^^ > m'jj^ 
above the lines and m^^ < m'jj^ below the lines. Also shown by the dashed lines are the mass 
contours for the tree- level value of mff±. The shaded region corresponds to the excluded region 
with < 0. 



^ 80 GeV [33, 34] limits us to the right of the m^± = 80 GeV contour. Therefore, 
requiring rrih^ > ttih^ while satisfying the experimental charged Higgs bounds restricts us 
to two regions: large A > 0.5, tua ^ ttiz for small tan /3 ^ 1 — 2, or small A < 0.5, tua ^ tuz 
for tan /3 ^ 2. Imposing a stronger bound on m^± from t bH^ searches at the Tevatron 
and the LHC [41-43] further narrows down the ruh^ > tuh^ region, resulting in a fine-tuned 
region to realize. 

On the other hand, rrih^ < tuh^ is much easier to realize in the NMSSM. In contrast 
to the MSSM, where being deep into the decoupling region rriA ^ 300 GeV is typically 
required to satisfy both the mass window and the cross section requirement (i.e. for to 
obtain SM-like couplings to the gauge bosons), in the NMSSM, with the mixture of the 
singlet and the possible suppressed couplings to 66, even a suppressed coupling to the gauge 
sector could be accommodated while satisfying the experimentally observed cross section 
range. Note that our discussions are based on tree level expression for the Higgs masses. 
While including loop corrections shifts all the masses, our statements are still qualitatively 
valid. 

Including the extra singlet in the spectrum gives three distinct cases, as sketched in 
Fig. 1, corresponding to either i^i, 7^2, or being SM-like: 

• Hi SM-like: rrih, < rriH^.ms, 

• H2 SM-hke: ms < rrih, < uih, or uih, < rrih, < ms, 

• SM-like: rriH^.ms < rrih^. 
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Table 1. Reduced Higgs couplings at tree level. The charged Higgs couplings of H^dLU^ and 
H~ULd^ are normalized to the SM top and bottom Yukawa couplings \^mt/v and V^m^/v, 
respectively. 

With the off-diagonal mixing in the mass matrix, the separation of these regions becomes 
less distinct while the above relations still approximately hold. 

2.3 Couplings 

The mass eigenstates i?i,2,3 are, in general, a mixture of /i^, Hy^ and S: 

Hi = Y.^H:Ha, fov i = 1,2,3, H^ = {K,H^,S), (2.13) 

a 

with ^^"^ being the 3x3 unitary matrix that rotates the Higgs bosons into the mass 
eigenstates. In particular, 1^^°"^ defines the fraction of hy^ Hy, and S in Hi with the 
unitarity relations: 

le^ri' + \^h;\' + \& = 1' \^Hf + i^Si' + \^Hf = 1- (2-14) 

Similarly, for the CP-odd Higgs bosons, the unitary rotation is Ai = J^a^A.'^^f^ where 
z = 1,2, and = (Ay, As). The fractions of Ay and As in the CP-odd mass eigenstates 
Ai,2 are given by l^^^p, i ^ 1,2, with = I^^P = 1 - l^fj' = 1 - 1^^ I'- 

In Table 1, we express the tree-level reduced couplings of the NMSSM Higgs mass 
eigenstates to various pairs of SM particles, which are the ratios of the NMSSM Higgs 
couplings to the corresponding SM values. The charged Higgs couplings of H^dLU^^ and 
H~ULd^^ are normalized to the SM top and bottom Yukawa couplings \/2mt/v and v^m^/i;, 
respectively. In the NMSSM, the HiZZ and HiWW couphngs are always modified in 
the same way at leading order. Therefore, we use VV to represent both WW and ZZ. 
The coupling of the CP-even Higgs bosons to the gauge boson sector VV is completely 
determined by the /i-^-fraction of Hi'. \^^.\'^^ which plays an important role in understanding 
the coupling and branching fraction behavior of the SM-like Higgs boson. Note that |^^. p < 
1, therefore, the HiVV couplings, as well as the Hi VV partial decay widths, are always 
suppressed compared to their SM values. However, the branching fractions of Hi VV 
could still be similar or even enhanced compared to their SM values, since Hi bb could 
be suppressed as well. 

The Higgs to 77 and Higgs to gg couplings are both loop-induced. The dominant 
contribution to the /i^77 coupling comes from the WW loop, with a sub-leading destructive 
contribution from the top loop. The hygg coupling, on the other hand, is dominated by the 
top-loop contribution. The i?i77 and Higg couplings are modified similarly in the NMSSM, 
based on the reduced couplings as listed in Table 1. 
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3 Parameter Scan and Constrained Regions 



We will focus our scan on the parameters that are most relevant to the Higgs sector, namely, 
parameters appearing in the Higgs potential, as well as the stop mass parameters, which 
could induce a relatively large loop correction to the Higgs mass. Since the impact of other 
SUSY sectors to the Higgs mass is typically small, we effectively decouple them by setting 
all other SUSY mass parameters to be 3 TeV and the other trilinear soft SUSY breaking 
parameters to be 0. Note that the sbottom and stau might modify the Higgs mass and 
certain couplings at loop level, which could have substantial effects in certain regions of 
parameter space. We defer a discussion of these regions to specific studies in the literature 
[44] and will only focus on the Higgs and stop sectors in the current study. 
In the MSSM, the relevant Higgs and stop sector parameters are 

mA, tan/3, /i, M^sQ, M^su, A, (3.1) 

as well as the Higgs vacuum expectation value v = 246 GeV. In the NMSSM, the tree level 
Higgs potential involves seven parameters: (A, a^, A^, A/^, '^s, tan/3) and v. After replacing 
Vs hy fj. = Xvg/V^ and replacing Ax by rriA as defined in Eq. (2.4), we are left with three 
new parameters compared to the MSSM case. We scan these parameters in the range of 



< 


mA 


< 200 GeV, 


1 < 


tan P 


< 10, 


100 GeV < 




< 1000 GeV 


0.01 < 


A, K 


< 1, 


1200 GeV < 


An 


< 200 GeV, 



100 GeV < MssQ, Mssu < 3000 GeV, 
-4000 GeV < At < 4000 GeV, (3.2) 

unless otherwise stated. The range of rriA is chosen to guarantee that Hy and Ay are light. 
The ranges of /x. A, and A,^ are chosen such that the CP-even and odd singlet masses are 
allowed to vary over a wide range. The stop sector mass and mixing parameters are chosen 
to cover both the minimal and maximal mixing scenarios. We restrict tan/3 to be in the 
range of 1 — 10 since regions with higher values of tan /3 do not contain a SM-like Higgs 
boson in the mass window of 124 — 128 GeV, as will be discussed in detail in Sees. 4 and 5. 

The scan is performed by utilizing NMSSMTools 3.2.1 [45-47] to calculate the Higgs 
and SUSY spectrum, Higgs couplings, decay widths, branching fractions, and various Higgs 
production cross sections. The full constraints imposed for the scan procedure include: 

• the latest LHC limits in various SM Higgs searches [48-54]; 

• bounds on MSSM Higgs search channels from LEP, the Tevatron, and the LHC [14, 
55]; 

• stop and sbottom masses to be heavier than 100 GeV. 
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We did not impose bounds that are not directly relevant to the Higgs sector, for example, 
other SUSY particle searches, flavor physics, and dark matter relic density. Those bounds 
typically involve SUSY parameters of other NMSSM sectors which we did not scan. Al- 
though some significant reduction of the allowed parameter space may occur with these 
additional constraints, we do not expect our conclusions to be changed. We generated 
a large Monte Carlo sample over the multi-dimensional parameter space and tested each 
parameter point against the experimental constraints. For the following presentation, the 
allowed points (or regions) in the plots are indicative of consistent theoretical solutions 
satisfying the experimental constraints, but are not meant to span the complete space of 
possible solutions. 

Given the discovery of a SM-like Higgs boson around 126 GeV, we study its implication 
by applying the following requirements step by step: 

Either Hi, or H2, or in the mass window of 124 - 128 GeV, (3.3) 

a X Biigg ^ Hi ^ ax Biigg ^ Hi ^ WW/ZZ) ^ , , 

/^^ ^ , ' ^ > 0.8, ^ ^ > 0.4. 3.4 

(cr X Br)sM {o- X Br)sM 

The cases delineated in Eqs. (3.3) and (3.3) determine the defining feature of the regions 
described in this paper and will, henceforth, be referred to as i7i-126, H2-I2Q and 7^3-126, 
respectively. 

Figure 3 shows the allowed mass regions versus m/^± for the CP-even Higgs bosons 
(left panels) and the CP-odd/charged Higgs bosons (right panels). The first and second 
row panels are for points that pass all the experimental constraints as itemized earlier, as 
well as Hi and 7^2, respectively, satisfying both the mass and cross section requirements as 
listed in Eqs. (3.3) and (3.4). The third row panels are for points that pass all experimental 
constraints as well as Hs satisfying the mass requirement as listed in Eq. (3.3). We have 
chosen to plot the physical Higgs masses against the charged Higgs mass m^± , rather than 
the conventional choice of tua as in the MSSM. Due to the relatively large loop corrections 
to the Higgs masses, the natural scale parameter choice in the NMSSM would be the loop 
corrected Ay mass ^A^^op' NMSSM equivalent of the MSSM ttia (the physical mass 
for the CP-odd MSSM Higgs). Ay, however, has to mix with As to provide masses for 
the two CP-odd mass eigenstates Ai and A2. The charged Higgs mass m^±, on the other 
hand, retains roughly the simple relationship with rriAi^^p^ described in Eq. (2.6), after loop 
corrections. Therefore, we choose the physical m^± as the scale parameter in Fig. 3. In 
this figure, we scanned in the range GeV < < 300 GeV, rather than GeV < rriA < 

200 GeV, to improve the coverage of the parameter region of our interest. 

For Hi being the SM-like Higgs in the mass window of 124 — 128 GeV (see the first 
row of Fig. 3), H2 is typically in the mass range of 125 — 300 GeV, while niHs ^ 200 GeV. 
The charged Higgs mass is in the approximate range 125 — 300 GeV. Charged Higgs bosons 
with mass less than 150 GeV are mostly ruled out by the direct search for H^ produced 
in top decays. The light CP-odd Higgs could be very light, a few GeV < rriAi ^ 300 GeV, 
while 771^2 ^ 200 GeV. When rriAi < mif^/2, the decay channel Hi AiAi opens, leading 
to very interesting phenomenology, as will be discussed in detail in Sec. 4. Note that the 
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Figure 3. The left panels show the allowed mass regions versus for the CP-even Hi (red), 

H2 (green), and Hs (blue). The right panels show the allowed mass regions versus m//± for the 
CP-odd Ai (magenta), and A2 (brown), as well as the charged Higgs (cyan). The first and 
second row panels contain the points that pass all the experimental constraints as well as having Hi 
and i^2, respectively, satisfying both the mass and cross section requirements as listed in Eqs. (3.3) 
and (3.4). The third row panels contain the points that pass all the experimental constraints as 
well as Eq. (3.3), but not Eq. (3.4). 
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boundary of the H2 and regions, as well as the boundary of the Ai and A2 regions show 
nice correlation with m^±. This is because the boundary is given by ttih^ as in Eq. (2.10) 
for the CP-even case, and by ttia for the CP-odd case, both of which scale with m^±. The 
singlet mixing with Hy and will push/pull the mass eigenstates away from ttih^ and 
m^, leaving a clear boundary. Given a 7^2, ^^3 pair (Ai, A2 pair), the one whose mass is 
closer to the H2-HS (^1-^42) boundary line is more Hy {Ay)-like. 

For H2 being the SM-like Higgs in the mass window of 124 — 128 GeV (see the second 
row of Fig. 3), a large fraction of the points contain Hi in the mass range of 60 — 124 GeV. 
There is also a significant set of points with tuhi < m/f2/2, which turns on the decay channel 
H2 HiHi^ as will be discussed in Sec. 5. tuhs is in the mass window of approximately 
200—350 GeV, and grows roughly linearly with m^±^ an indication of i^a being mostly 
iJ-y-like. The points with ttih^, below ^180 GeV are removed by a combination of the 
collider constraints and the cross section requirement of Eq. (3.4). This is very different 
from the i7i-126 case, in which could be singlet dominant with mass as large as 1 TeV 
or higher. For the light CP-odd Higgs Ai, it falls into two regions: one region with 60 GeV 
^ ^ 300 GeV {m^± > 200 GeV), with little dependence on m^± (for Ai being mostly 
^5); another region with ttiai ^150 GeV {ttih± ^150 GeV), which grows linearly with 
m^± (for Ai being A^-like). A2 typically has a mass of 200 GeV or higher, which also falls 
into two regions accordingly. 

For i^a being the SM-like Higgs in the mass window of 124 — 128 GeV (see the third row 
of Fig. 3), both the singlet and iJ-^-dominant Higgs bosons need to be lighter than about 126 
GeV. Given the tight experimental constraints on the light Higgs searches, as well as the fine- 
tuning between the mass parameters, this region turns out be to highly restrictive. While we 
can realize regions with ttih^ in the desired mass window, it is extremely difficult to satisfy 
the cross section requirement of Eq. (3.4). Panels in the third row of Fig. 3 show points with 
in the mass window of 124 - 128 GeV. However, a x Br(^^ ^ H3 ^ -f-f, WW/ZZ)/SM 
is less than 0.4 in general, and it is therefore hard to accommodate the observed Higgs 
signal as Hs in the NMSSM. 

In what follows, we will discuss the i7i-126 and H2-I2Q cases in detail, exploring the 
relevant parameter space for each region, the composition of the 126 GeV SM-like Higgs 
and the other light NMSSM Higgs bosons, possible enhancement or suppression of various 
search channels as well as correlations between them. 

4 Hi as the SM-like Higgs Boson 
4.1 Parameter Regions 

For Hi to have SM-like cross sections for gg Hi 77, WW/ZZ within the experimen- 
tally observed ranges. Hi needs to be either dominantly hy or have a considerable singlet 
fraction with a suppressed Hi bb partial width. Hy and S dominant states are typi- 
cally heavier such that, usually, the lightest CP-even Higgs state is mostly hy. This case is 
seldom realized in the MSSM low-m^ region (m^ < 2m^), since the light CP-even Higgs 
boson typically has suppressed couplings to WW and ZZ in this region. In the NMSSM, 
the tree-level diagonal mass term for hy is m^^ = cos^ 2/3 + ^{Xv)'^ sin^ 2/3. Large A and 
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Figure 4. The dependence of ttihi on the fohowing NMSSM parameters in the i^i-126 case: A, 
hz, tan /3, and /i. Grey points are those that pass the experimental constraints, pale-pink points are 
those with Hi in the mass window 124 GeV < tuhi < 128 GeV and green points are those with the 
cross section requirements further imposed. Black points are those that remain perturbative up to 
the Planck scale. 



small tan^ are preferred to push up the mass of hy into the desired mass window. For 
small tan/3, even for small m^, typically m^^ < Toj^^, resulting in the lighter MSSM-like 
CP-even Higgs being SM-like in the low-m^ region. In addition, mixture with the singlet 
in the NMSSM which produces, in particular, a suppressed Hi bb partial decay width, 
could lead to a SM-like 77 and WW/ZZ branching fraction for Hi as well. The push-down 
effect in mass eigenvalues from the singlet mixing also helps to realize the mostly hy state 
being Hi. 

To show the effect of the narrowing down of the parameter regions due to the mass and 
cross section requirements. Fig. 4 presents the dependence of tuhi on A, ^c, tan/3, and /x, 
with gray dots for all points satisfying the experimental constraints, pale-pink points which 
pass the mass window requirement of Eq. (3.3), and green points, that almost overlap the 
pale-pink points, which pass both the mass and cross section requirements of Eqs. (3.3) and 
(3.4). 

After requiring Hi to fall into the mass region of 124 — 128 GeV, we are restricted 
to the parameter region of A ^ 0.55, k ^ 0.3 (with a small number of points down to 
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(GeV) 


Hi-126 
perturb. 


1-3.5 
1.5 -2.5 
1-3.5 


0-200 
150-200 
100-200 


100-500 
100-150 
100-200 


> 0.55 
0.55-0.65 

> 0.55 


^ 0-3 
0.3-0.5 
>0.5 


-650-300 
-30-230 
-150-150 


-1200-200 
-150-100 
-50-30 



Table 2. NMSSM Parameter regions for the i^i-126 case. 

0.1), 1 < tan/3 < 3.5, /a < 500 GeV, -1200 GeV < A^ < 200 GeV with no restriction 
on rriA which is allowed to be in the entire region of — 200 GeV (the corresponding 
region for Ax is approximately —650 GeV to 300 GeV). The stop mass parameters M^sq^ 
M^su and At are unrestricted as well. Further imposing the cross section requirement 
for gg ^ Hi ^ 77, WW and ZZ does not narrow down the allowed regions for these 
parameters further. 

Also shown as the black points in Fig. 4 are the parameter points where A and k remain 
perturbative up to the Grand Unified Theory (GUT) scale. They occupy a small region of 
0.5 < A < 0.65, 0.3 < < 0.5, tan/3 - 2, 100 < ^ < 150 GeV, -150 < A^ < 100 GeV, and 
150 <mA< 200 GeV (-30 < Ax < 230 GeV). While M^sq and M^su are unconstrained, 
\At\ is restricted to be > 1200 GeV. These parameter regions are summarized in Table 2. 

We have noted earlier that the light CP-odd Higgs Ai could be very light. When it 
falls below half of the Hi mass. Hi AiAi opens up, which could dominate the Hi decay 
width, compared to the usual case in which decay to bb dominates. Therefore, we further 
separate the i7i-126 case into three regions: 

• Hi Region lA: m^i > and > 0.7: green points in Figs. 5-10. 

• Hi Region IB: m^i > and |C^^J^ < 0.7: red points in Figs. 5-10. 

• Hi Region II: m^i < mif^/2: magenta points in Figs. 5-10. 

To identify the NMSSM parameter regions that give a SM-like Hi in the mass window 
of 124 — 128 GeV, in Fig. 5, we show the viable regions in various combinations of NMSSM 
parameters. Grey points are those that pass the experimental constraints, pale-pink points 
are those with Hi in the mass window 124 GeV < tuhi < 128 GeV, green and red points 
are for Hi Region I with tuai above the Hi AiAi threshold, and magenta points are 
for Hi Region II with low m^^. Again, the black points are those where A and n remain 
perturbative up to the GUT scale. 

The first two panels show the (a) A versus m^, and (b) n versus ttia regions. For small 
values of m^, A has to be around 0.6—0.7, since too large a value of A is ruled out by the 
charged Higgs mass bounds, while too small a value of A results in ttihi being less than 124 
GeV. For larger values of m^, the A range is enlarged to 0.55 ^ A < 1. on the other 
hand, has to be ^ 1 for small m^, while smaller n is allowed for larger m^. 

Panel (c) of Fig. 5 shows the viable region in the A versus n plane. Given A > 0.55 
and Kj > 0.3, the renormalization group running of A, and as well as of the Yukawa 
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Figure 5. Viable NMSSM parameter regions in the Hi-126 case: (a) A versus m^, (b) hi versus m^, 
(c) A versus hz, (d) tan/3 versus rriAi^^p, (e) Ax versus /i, and (f) versus /i. Grey points are those 
that pass the experimental constraints, pale-pink points are those with Hi in the mass window 124 
GeV < rriHi < 128 GeV. Green points are for Hi Region lA: m^i > mnj^ and |^^"J^ > 0.7. Red 
points are for Hi Region lA: m^i > ttihi/'^ and I^^^^J^ < 0.7. Magenta points are for Hi Region 
II: iriAi < mHif^- The black points are those where A and remain perturbative up to the GUT 
scale. 



couplings yt^b and gauge couplings might reach the Landau pole before Mqut- As noted in 
Ref. [6], a larger A allows a highly natural light Higgs boson. For all the points that pass 



the mass and cross section requirements, only a small region of the X-hz plane, as shown by 
the black points in panel (c) of Fig. 5, remain perturbative up to the GUT scale around 
10^^ GeV. For larger values of ^c, it reaches the Landau pole before the other couplings. 
While the running of A is much slower, it has a large impact on the running of the gauge 
couplings and Yukawa couplings. Increasing the value of A would accelerate the running 
of the top Yukawa coupling. However, for all the viable points that pass both the mass 
and cross section requirements, the scale at which at least one of the couplings becomes 
non-perturbative is typically larger than 10^ GeV, much higher than the electroweak scale. 
Since adding new multiplets or other new physics could affect the running of the couplings 
and delay the Landau pole scale, in our study, we relax the perturbativity constraint and 
only place a loose upper bound of A, < 1. All of our parameter points remain perturbative 
up to at least the scale of 10^ GeV. 

Panel (d) of Fig. 5 shows the viable region in the tan/^-m^^^^^ plane, where we have 
plotted rriAi^^p for better comparison with the MSSM.^ Unlike in the MSSM case, where 
constraints from collider direct Higgs searches and the light CP-even Higgs being SM-like 
require the parameters to be in the decoupling region of rriA ^ 300 GeV [16], in the NMSSM, 
by contrast, with the push-down effect from the singlet mixing and the extra contribution 
from ^(Ai;)^ sin^ 2/3 to the tree level mass squared for /i^. Hi could be the SM-like Higgs 
in the low-m^ region: rriAi^^p ^ 140 GeV (while the tree-level rriA could be as low as a few 
GeV). The range of tan/3 is smaller compared to that of the MSSM:^ 1 < tan/3 < 3, since 
smaller tan /3 is preferred for providing a sizable contribution from the A-term to the tree 
level Higgs mass m^^ . 

Panel (e) of Fig. 5 shows a clear correlation between Ax and /x. This is because a larger 
value of fi is needed to cancel the negative contribution from Ax to keep m\ > 0, as given 
in Eq. (2.4). Panel (f) of Fig. 5 shows a weaker correlation between A,^ and fi. While larger 
fi is typically preferred for a larger negative A^, can not be too large since otherwise at 
least one of the CP-even Higgs masses squared becomes negative. 

The magenta points in Fig. 5 are in Hi Region II: ruAi < mif^/2. It maps out the region 
of smah \A^\, \Ax\ and /i: -50 < A^ < 30 GeV, -150 < Ax < 150 GeV, 100 < /x < 200 
GeV. rriA is restricted to be in the range of 100 — 200 GeV {jtiai^^^ ^150 GeV), and in 
the range of 0.5 — 1. Ranges for A and tan/3, however, are not narrowed compared to the 
generic Hi Region I, as shown as the green points in Fig. 5. 

Unlike in the MSSM case, where the mass parameters for the stop sector, M^sQ and 
Mssui are correlated with the stop left-right mixing At to be close to the m^^^ scenario, 
l^ltl ^ ^JQM^SQM^SU 1 there is no obvious correlation between M^sQ^ M^sui and At in the 
NMSSM. All the ranges are allowed for these parameters. This is because in the MSSM, 
we need large loop corrections to the Higgs mass from the stop sector to push it to the 
124 — 128 GeV mass window, which requires either large stop masses around 5 — 10 TeV 
or large stop mixing. In the NMSSM, such a lift to the Higgs mass could be achieved by 
the (Xv)'^ contribution to the Higgs mass at tree level, resulting in a less constrained stop 

^niA in the MSSM is the physical mass for the CP-odd Higgs with loop corrections already included. 
^tan/3 < 3 is excluded by the LEP Higgs searches in the MSSM [60]. 
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Figure 6. The normalized a x Br/SM for (a) gg ^ Hi ^ 77 and (b) gg ^ Hi ^ WW/ZZ as 
a function of tuhi in the i^i-126 case. The current experimental constraints from the SM Higgs 
searches of the 77, WW and ZZ channels are also imposed. Color coding is the same as for Fig. 5. 



sector. The mass for the hghtest stop can be as hght as 100 — 200 GeV, with shghtly larger 
mass splitting Am^ > 200 — 300 GeV anticipated for small m^^ . However, once m^^ > 400 
GeV, a degenerate stop mass spectrum can also be accommodated. 

4.2 Production Cross Sections and Decay Branching Fractions of the SM-hke 

Hi 

In Fig. 6, we show the cross sections ratios of NMSSM model to the Standard Model 
(cr X Br/SM) for gg ^ Hi ^ 77 in panel (a) and gg ^ Hi ^ WW/ZZ in panel (b) as 
a function of ttihi- The current 95% C.L. experimental exclusion limits for the SM Higgs 
searches in the 77, WW and ZZ channels are also imposed. While the 77 limit imposes 
strong constraints in the low rriHi region, for ttihi in the mass window of 124 — 128 GeV, 
the WW and ZZ cross section bounds are more important and rule out points with large 
a X Br. For the 77 channel, a x Br/SM mainly varies in the range of 0.8 — 1.75, where 
the lower limit comes from our requirement of the signal region, as indicated by the current 
Higgs signal at both the ATLAS and CMS experiments [4, 5, 51]. Notice that for a few 
points, a x Br/SM as large as 2 can be reached. For the WW/ZZ channel, a x Br/SM 
varies mostly between the range of 0.7 — 1.6, with a few points that could reach a value of 
0.5 or even smaller. 

In the NMSSM, both the production cross section and decay branching fractions could 
deviate from their SM values. In the mass window of 124 — 128 GeV, cr(gg Hi) /asm 
typically varies between 0.6 — 1.4, although a suppression as small as 0.2 or an enhancement 
as large as 1.7 are also possible. For the decay branching fraction. Hi WW/ZZ (77) is 
typically approximately 0.6—1.5 (0.6—2) of the SM value. There are a few points with very 
large enhancement factors, approximately 3—4 (5—6) for WW/ZZ (77), which are needed 
to compensate the associated suppression from the gluon fusion production. 

Hi 77 and Hi WW/ZZ are highly correlated, as in the case of the MSSM sce- 
nario. This is because the loop generated i?i77 coupling receives its dominant contribution 
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Figure 7. The normalized a x Br/SM for (a) 77 versus WW/ZZ channel and (b) the normalized 
Br/BrsM for WW/ZZ versus bb in the i^i-126 case. Color coding is the same as for Fig. 5. 




Figure 8. The normalized a x Br/SM for gg ^ Hi ^ 77 with (a) ttia dependence and (b) ji 
dependence in the i^i-126 case. Color coding is the same as for Fig. 5. 



from the VF-loop and is therefore controlled by the same HiWW coupling. Such correlation 
is shown in Fig. 7, panel (a) for 77 versus VV. In the Hi-126 case, most of the points fall 
into the region of 



a X Br(^^ ^ Hi^ 77) /(^ x Br)sM ^ Br(gi 77)/BrsM 
a X Br(^^ ^Hi^ VV)/{a x Br)sM ~ Bt{Hi ^ VV)/Btsm 



1.1 . (4.1) 



However, there are a few scattered points with larger 77 : VV ratios. These points have an 
enhanced Hi 77 partial width due to light stop contributions. 

Unlike the correlation shown in the 77 versus VV channel, the correlation between 
the bb and VV channels exhibit interesting feature, as shown in Fig. 7, panel (b), for 
bb) versus b^(^i ~^ VV). For Hi Region I, bb and WW are anti-correlated 
so that the VV channel is enhanced compared to the SM value only when the bb channel 
gets relatively suppressed. This is as expected since bb and VV are the two dominant Hi 
decay channels for rriHi > 2m^^ . 
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For Hi Region II with low m^^, however, no such correlation is observed. While 
Hi bb could be much suppressed compared to its SM value, Br(i7i VV) /Btsm varies 
in the range of 0.5 — 1.5, almost independently of Br(i7i 66)/BrsM- The opening of 
the Hi AiAi channel in this mass window replaces Hi bb to keep Hi VV in the 
desired range to satisfy the cross section requirement. 

The TT and bb channels have also been searched for at the LHC, which indicate a weak 
SM Higgs signal of approximately 1 — 2a [49, 50, 53, 54]. For the rr channel, while the 
dominant contribution comes from the vector boson fusion (VBF) production, gg ^ H ^ 
TT could be separated with a dedicated search [49, 53]. H ^ bb has been studied for both 
VH and UH production, with better limits coming from VH associated production [50, 54]. 
In the NMSSM, since it is the same down-type Higgs H^ that couples to both the bottom 
quark and the tauon, Hibb and Hitt receive the same corrections (up to small difference 
in the radiative corrections that are non-universal for bottom and tau). Therefore, the 
bb and rr channels are highly correlated: Br(i7i 66)/BrsM ^ Br(i7i TT)/BrsM- 
For VBF/VH ^ Hi ^ TT/bb, a x Br/SM is < 1.1. While for gg ^ Hi ^ t^t' , 
an enhancement as large as 1.5 of the SM value is possible, which is again from stop loop 
corrections to gg Hi. UHi with Hi bb receives little enhancement, a x Br/SM < 1.05. 

Fig. 8 shows the parameter dependence of a x Br/SM for gg ^ Hi ^ 77 for ruA [panel 

(a) ] and /i [panel (b)]. Larger values for a x Br/SM is achieved for larger values of uia and 
smaller values of fi. If a significant enhancement of gg ^ Hi ^ 77 is observed in future 
experiments, tua and /i (as well as Ax) would be restricted to a narrower region. 

4.3 Wave Function Overlap 

The deviation of the production and decay of Hi can be traced back to the hy^ Hy and 
S fractions in i^i, which is given by the wave function overlap IC^^'J^, IC^J'P and |<^^J^, 
as defined in Eq. (2.13). Fig. 9 shows IC^.P versus \^^.\'^ for Hi [panel (a)], H2 [panel 

(b) ] and H^ [panel (c)]. Since + + \^fjf = 1, the distance between the cross 
diagonal line and the points indicates the value of iCjjP- Fo^* the generic Hi Region I A 
(green points), l^^""^ p + |^^^ P ^ 1; the iJ-^-fraction in Hi is almost 0. Typically, about 70% 
or more of Hi is /i^, which couples exactly like the SM Higgs, while the singlet component 
varies between to approximately 30%. For 7^2, it could be either i7-i;-dominant for those 
points with IC^^I^ ^ 0, or a mixture of Hy and S for points with larger |C^2l^' ^3 mostly 
singlet-dominant, or with a small mixture of hy for points close to the cross-diagonal line. 
It could also have a significant Hy — S mixture for points with smaller iC^gP- 

For the Hi Region IB (red points), \^^^\'^ < 0.7, Hi typically has a sizable fraction of 
Hy and S. H2 is also a mixture of hy^ Hy and 5, with H^ being mostly a Hy-S mixture. 

For Hi in the low rriAi region (magenta points of Hi Region II), Hi and H2 are mostly 
di hy — Hy mixture, with Hi being more hy-like, and H2 being more i^-^-like. This region 
and Region IB share the property that they typically depend on a suppressed Hibb coupling 
proportional to — ^^"^ tan/3. The S fractions of Region II vary between to 25% for 
both Hi and i72, while it is the dominant component of H^. 

Fig. 10 shows the fraction of MSSM CP-odd Higgs Ay in the light CP-odd Higgs Ai: 
I^^^P as a function of A,^ [panel (a)] and [panel (b)]. The more negative A,^ becomes. 
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Figure 10. A^-fraction in the light CP-odd Higgs Ai\ l^^^'p = 1 - = 1 - in the 

i^i-126 case. Color coding is the same as for Fig. 5. 



the larger the diagonal mass term m^^for the singlet Ag becomes, which results in Ai 
becoming more and more A^-like. The ji dependence also shows a trend of large ji leading 
to Ai being more A^-like, mainly due to the —3k/X jiA,^ contribution to the m^^ mass 
term. The A^-fraction in Ai, as well as the A-^-fraction in A2^ is simply 1 — iC^^p. For the 
points that satisfy the perturbativity requirement (black points), Ai is mostly singlet like. 
For regions with small tuai < rriHi/^ (magenta points), a significant fraction of Ai, 40% to 
80%, is singlet. While for generic Hi Region I, Ai could be either A-^-like (small m^, large 
negative A,^ and Ax^ large /x, large hc ) or A^-like, depending on the NMSSM parameters. 
Note that while we are focusing on the low ruA region, which controls the mass scale for the 
MSSM-type CP-odd Higgs, the mass parameter for the CP-odd singlet Higgs could vary in 
a large range given the scanning parameter region. As a result, rriAi is below 300 GeV, 
while the mass for the heavy CP-odd Higgs, mA2i could be as large as 1 TeV or higher. 
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rriA 
(GeV) 


(GeV) 


A 




Ax 
(GeV) 


(GeV) 


H2-I26 
perturb. 


^1 
1-3.25 
1.5-2.5 
1.25-2.5 


0-200 
100-200 
170-200 
125-200 


100-300 
100-200 
100-130 
100-150 


0-0.75 
0.4-0.75 

0.5-0.7 
0.5-0.75 


0-1 
> 0.05 
0.05-0.6 
>0.3 


-600-300 
-300-300 
0-300 
0-200 


-1200-50 
-1200-50 
-300-50 
-500—250 



Table 3. NMSSM Parameter region for the H2-I26 case. 

5 H2 as the SM-like Higgs Boson 
5.1 Parameter Regions 

In the hmit where the mixings between hy, Hy and S are smah, there are two cases that 
give rise to H2 being SM-hke: Hi being singlet hke and Hs being mostly Hy] or Hi 
being mostly Hy-like and Hs being mostly singlet. Including the loop corrections as well 
as mixture between hy^ Hy and 5, the separation between these two cases becomes less 
distinct. The former case is similar to i7i-126, except that the singlet is now the lightest 
state. The latter case is similar to the MSSM non-decoupling region, which requires a high 
level of fine-tuning to satisfy the experimental constraints, as well as mass and cross section 
requirements. As a result, in the NMSSM, while there are points with a relatively large 
fraction in Hi, there is always a sizable 5- fraction in Hi as well. 

Unlike the Hi-126 case, where imposing the mass window on rriHi already greatly 
narrows down the parameter regions while the cross section requirement usually does not 
provide further restriction, imposing the mass window in the H2-I26 case (124 GeV < 
m/f2 < 128 GeV) does not greatly reduce the parameter space beyond the already restricted 
space from satisfying the experimental constraints. Requiring H2 to have a SM-like gg 
77, WW/ZZ rate, however, does further reduce the parameter space to be in the small 
tan/3, small /x, medium to large A, and small \ Ax\ region, as summarized in Table. 3. Note 
that compared to the i7i-126 case, where tua could be very small, in the H2-I2Q case, ttia is 
typically larger than about 100 GeV. We note, however, that rriAi^^p in both cases is greater 
than approximately 150 GeV. In the Hi-126 case, the SM-like Higgs is pushed down and 
requires a larger stop-loop correction while in the present H2-I26 case, the SM-like Higgs 
is pushed up and, as a result, requires less of a contribution from the stop sector. The stop 
mass parameters M^sQ^ M^su and At are therefore less restricted in the H2-I2Q case. 

Also shown in Table 3 is the region where A and Hi remain perturbative until the GUT 
scale. Unlike the i7i-126 case in which l^^l is restricted to be > 1200 GeV, for the H2-I2Q 
case, \At\ is typically unrestricted. 

While the light CP-odd Higgs Ai is almost always heavier than 771/^2/2, the lightest 
CP-even Higgs Hi could be lighter than tuh^/'^ such that the H2 HiHi decay opens up. 
Although H2 is typically /i^-like, it could obtain a relatively large S^-fraction to suppress 
the otherwise dominant decay mode H2 bb. Therefore, we separate the H2-I26 case into 
three regions: 

• H2 Region lA: rriHi > ^if2/2 and l^^"" p > 0.5: green points in Figs. 11-14. 
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• H2 Region IB: rriHi > ^7^2/2 and < 0.5: red points in Figs. 11-14. 

• H2 Region II: ttihi < mif2/2: magenta points in Figs. 11-14. 

In Fig. 11, we show the viable regions in various combinations of the NMSSM param- 
eters. The first two panels show the (a) A versus m^, and (b) k, versus rriA regions. Small 
ruA ^ 100 GeV is not favored since the cross sections for gg ^ H2 ^ ^^,WW/ZZ are 
suppressed. A is typically in the range of 0.4 — 0.75. Smaller A is not allowed due to the 
suppressed cross sections, while larger values of A are not allowed due to charged Higgs 
mass bounds, k varies over the whole range of — 1, with larger values of uia preferred 
for smaller k. 

Panel (c) of Fig. 11 shows the viable region in the X-hc plane. Regions with small A 
satisfy the mass window but fail the cross section requirement. Also, shown in black, are 
those points that remain perturbative until the Planck scale, which spans a range of A 
between 0.5 to 0.7 and Hi between 0.1 and 0.5. Panel (d) of Fig. 11 shows the viable region 
in the tan/3-m^^^^^ plane, tan/3 falls into a range of 1.5 — 3.25, while rriAi^^p varies between 
160 - 240 GeV (m^ varies between 100 - 200 GeV). 

Panel (e) of Fig. 11 shows a weak correlation between Ax and /i. Regions of Ax < —300 
GeV fail the cross section requirement. There is also a correlation between A,^ and /x, as 
shown in panel (f) of Fig. 11. This is because in the 7^2-126 case, most Hi are singlet- 
like. The CP-even singlet mass needs to be smaller than ruh^ and is typically controlled 
by the cancellation between a positive parameter and a negative A,^ term, as shown in 
Eq. (2.11). This correlation can be seen more clearly in H2 Region II (magenta points) 
where finer cancellation is enforced. 

The magenta points in Fig. 11 are for H2 Region II: rriHi < 771/^2/2. They span the 
region of small \A,^\^ \Ax\ and /x, intermediate ^c, tua ^ 200 GeV, and 1.5 < tan/3 < 2, as 
summarized in Table 3. 

5.2 Production Cross Sections and Decay Branching Fractions for the SM-hke 

The ranges of a x Br/SM for gg ^ H2 ^ 77, shown in Fig. 12, is slightly large than that 
of the i7i-126 case. An enhancement as large as a factor of 2 can be achieved in the present 
case. For gg ^ H2 ^ WW/ZZ, the a x Br/SM is typically in the range of 0.4 - 1.6, and 
bounded above by the current experimental searches in the WW/ZZ channels. Note that a 
relatively strong suppression of about 0.4 could be accommodated more easily than in the 
Hi-126 case. 

H2 77 and H2 WW/ZZ are also highly correlated, as shown in Fig. 13, panel 
(a) for 77 versus WW/ZZ. There are several branches, corresponding to H2 Region I A 
and IB as categorized in Sec. 5.1. For Region lA (green points) with H2 being mostly hy- 
dominant, Br(^iff^^"^)/^^ ^ ^-^ lower branch of green points. However, there is 

another branch with a higher value of B^^lff ^^w^) / CBr )^m ^ ^' '^'^^^^ points typically have 
an enhanced H2 77 compared to the SM value due to the light stop contributions. For 
Region IB (red points) with H2 being a mixture of hy, Hy and S, ^^hI^^^w^^ ^ 
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(e) (f) 

Figure 11. Viable NMSSM parameter regions in the H2-I26 case: (a) A versus m^, (b) versus 
m^, (c) A versus hz, (d) tan/3 versus rriAi^^p, (e) Ax versus /i, and (f) versus fi. Grey points 
are those that pass the experimental constraints, pale-pink points are those with H2 in the mass 
window 124 GeV < < 128 GeV. Green points are for H2 Region lA: tuhi > mH2/'^ with 

I^^Hsl^ > 0.5. Red points are for H2 Region IB: rriHi > ^i^s/^ with l^^^l^ < 0.5. Magenta points 
are for H2 Region II: rriHi < ttih^/'^' The black points are those where A and n remain perturbative 
up to the GUT scale. 

In Fig. 13, panel (b), we show the correlation between the hh and VV channel: Br(i72 
66)/BrsM versus Br(i72 VV) /^isM- While most regions exhibit an anti-correlation as 
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Figure 12. The normalized a x Br/SM for (a) gg ^ H2 ^ 77 and (b) gg ^ H2 ^ WW/ZZ as 
a function of the H2-I26 case. The current experimental constraints from the SM Higgs 

searches of the 77, WW and ZZ channels are also imposed. Color coding is the same as for Fig. 11. 




Figure 13. The normalized a x Br/SM for (a) 77 versus WW/ZZ channel, and the normalized 
Br/BrsM for (b) WW/ZZ versus bb in the case. Color coding is the same as for Fig. 11. 



expected, in H2 Region II (magenta points) with tuhi < m/f2/2, the branching fraction for 
the VV channel is almost independent of the bb channel. This is, similar to the magenta 
region in the Hi-126 case, due to an opening up of the decay channel H2 i^ii^i, which 
compensates for the suppression of the bb channel while keeping the total decay width of 
H2 close to the SM value. 

The bb and rr channels also exhibit a similar correlation behavior as in the Hi-126 
case: Br(i72 bb)/BTsM ^ Br(i72 tt)/BtsM' For VBF and VH with H2 rr.bb, 
a X Br/SM is in the range of 0.4 — 1.1 for Region lA and is much suppressed in Region 
IB and is < 0.4. For gg ^ H2 ^ r^r~ ^ most of the H2 Region lA falls into the range of 
0.4 — 1.4, although an enhancement as large as 2 is possible. For Region IB, this channel 
is almost always suppressed with a x Br < 0.8(cr x Br)sM- The process ttH2 with H2 bb 
receives little enhancement, with a x Br/SM < 1.06 for Region I A and a x Br/SM ^0.7 
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I^Hll I^H2l I^HS' 

(a) (b) (c) 

Figure 14. |^^. P versus l^^"". P for JYi (a), H2 (b) and (c) in the i^2-126 case. Color coding 
is the same as for Fig. 11. 



for Region IB. 

5.3 Wave Function Overlap 

Fig. 14 shows IC^.P versus \^^.\'^ for (a) Hi, (b) H2 and (c) H^. For H2 Region II (magenta 
points), Hi is mostly singlet, H2 is mostly hy and Hs is mostly 

For H2 region lA (green points) with > 0.5, while H2 is mostly hy-like by 

definition, its i7-y-fraction is almost always small. In contrast, while Hi is dominated by 5, 
it could have a relatively large i7^-fraction. H^ is typically a mixture of S and Hy, with the 
i7-y-fraction always being sizable: iC^g P ^ 0.4. The /i-i;-fraction in H^ is almost negligible. 

For H2 region IB (red points) with < 0.5, the singlet fraction in H2 could 

be significant, sometime even as large as 0.8. While the /i-^-fraction in H2 decreases, it 
increases accordingly in Hi: \^^^\'^ > 0.5. This opens up the possibility of Hi with sizable 
HiWW/HiZZ couplings that we will discuss in the next section. Both Hi and H2 could 
have a fraction of Hy as large as 0.3— 0.4. i^a, on the other hand, is mostly a mixture of 
Hy and 5, with the /i-^-fraction being negligible. 

The compositions of Ai and A2 are similar to that of the Hi-12Q case. Larger negative 
values of A,^ lead to a large fraction of Ai being Ay. However, for A,^ ^ 0, Ai could be 
mostly Ag. 

6 LHC Phenomenology for the Non-SM-like Higgs Bosons 

In the previous sections, we have presented two very interesting scenarios in the low-m^ 
region. The SM-like Higgs boson could be the lightest scalar particle (i7i-126) while the 
next lightest one is an admixture of its MSSM partner and the singlet state. The alternative 
is that the SM-hke Higgs boson is the second hghtest {H2-I26) while the hghtest scalar is a 
Hy-S-hy mixture. While the collider phenomenology of the SM-like Higgs boson has been 
shown earlier, it would be interesting to identify the signal features of the other low-mass 
Higgs bosons. 

6.1 Hi as the SM-like Higgs Boson 

In Fig. 15(a), we show the dominant production cross sections of gg fusion and VBF for 
Hi (red and pink points), H2 (green and light green points) and H^ (blue and light blue 
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Figure 15. Cross sections at the 14 TeV LHC in the Hi-126 case for (a) i^i,2,3 production via 
gg fusion (VBF) denoted by red (pink), green (light green), blue (light blue) points, respectively, 
and for (b) ^41^2 production via gg fusion {bb fusion) denoted by purple (light purple), brown (light 
brown) points, respectively. The yellow lines indicate the cross sections with SM couplings. 



points), respectively, satisfying all the constraints for the Hi-12Q case at the 14 TeV LHC. 
The yellow lines indicate the corresponding cross sections with SM couplings. When the /in- 
fraction is sizable, the production cross sections for H2 could be similar to the SM-like rate. 
The cross sections could also be suppressed by two orders of magnitudes if the S^-fraction 
is large, as for the Hs case. The VBF process can be more significantly suppressed than 
that of gg fusion. The production cross section for the CP-odd states ^41^2 from gg fusion 
via triangle loop diagrams is shown in Fig. 15(b). The rate can be similar to that of the 
SM-like Higgs boson and the spread of the cross section over the parameter scan is roughly 
about an order of magnitude, less pronounced than those for the CP-even cases. Although 
about an order magnitude lower, the production cross section from bb annihilation can be 
significantly larger than that of the SM value, due to the tan/3 enhancement. 

In Fig. 16, we further show the relevant branching fractions for i?i,2,3, ^1,2 and 
to the SM particles (and 7^2,3 for the case of Hi-126. The yellow lines indicate 

the corresponding branching fraction values with SM couplings. The non SM-like Higgs 
bosons typically have suppressed decay branching fractions to the regular SM channels, in 
particular for H^^ due to the opening up of new decay channels to lighter Higgs bosons 
pairs. The experimental searches for those new Higgs bosons at the LHC will continue to 
cover a broad parameter region. Ai decays dominantly to 66, with about 10% to t+t~, as 
shown in Fig. 16(e). The phenomenological consequences of this decay have been studied 
in the literature [32, 56-58], emphasizing the h A\A\ 4t, 46, 2t26 modes and we will 
not discuss them further here. One of the most striking results for the CP-odd Higgs decay, 
perhaps, is the potentially very large enhancement for the branching fraction A\ 77, as 
seen in Fig. 16(e). This is partly because of the reduced F^^t caused by the suppression 
of the A^b coupling, and partly because of the enhanced T{^A\ 77) due to the loop 
contributions from the light charginos and charged Higgs bosons, and from the top and 
stop. In the pure singlet limit, the dominant viable decay channel is A\ 77 induced 
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Figure 16. Decay branching fractions for i^i,2,3, ^1,2 and to the SM particles (and i^2,3 hh) 
in the case of Hi-126. The yehow hnes indicate the corresponding values with the SM couplings. 

by the chargino loop and charged Higgs loop from their non-suppressed couplings with the 
singlet. However, the chargino in our case is always much lighter than the charged Higgs, 
granting non-zero coupling. The total width could be as low as around 10~^ GeV. 

This may lead to interesting scenarios with a proper LSP that produces a greatly suppressed 
low-end 77 continuum for an indirect dark matter search [59] such as Fermi-LAT. 

Another interesting feature is that the CP-even heavy Higgs bosons could decay to a 
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Figure 17. Decay branching fractions for i^i,2,3; 
Hi-126. 



A\^2 and to Higgs bosons in the case of 



pair of stops when kinematically accessible. It is important to note that a heavier Higgs 
boson could decay to a pair of lighter Higgs bosons at a substantial rate and sometimes 
dominantly, as long as kinematically accessible. As shown in Fig. 17 for the case of i7i-126, 
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Figure 18. Cross sections at the 14 TeV LHC in the H2-I26 case. The color codes and the legends 
are the same as in Fig. 15. 

we see that 



Hi^AiAi, ZAi, (6.1) 

H2^AiAi, ZAi, HiHi, (6.2) 

Hs^AiAi, HiHi, ZAi, W^H^, A1A2, H2H2, H^H', (6.3) 

H^^W^Ai, W^H2, W^Hi, (6.4) 

Ai^ZHi, (6.5) 

A2^AiHi, A1H2, W^H^, ZHi, ZH2, ZHs, A^H^, (6.6) 



roughly according to the sizes of the branching fractions at the low values of the mass. 
The relative branching fractions depend on phase space factors and the couplings dictated 
by the MSSM and singlet components. Consequently, the striking signals will be multiple 
heavy quarks, such as 46, 4t and 262t, and will likely include t^t~ as well. While the final 
state with a or Z may be a good channel from the event identification view point, the 
final states with multiple heavy quarks may be rather challenging to separate out from the 
large SM backgrounds. 

6.2 H2 as the SM-like Higgs Boson 

Similar results for the Higgs production and decay channels are shown in Figs. 18—20, 
respectively, at the 14 TeV LHC for the H2-I26 case. It is interesting to note that Hi is 
non-SM-like, and lighter than 7^2, yet it could have as large a production cross section as 
H2. Although the branching fractions to WW, ZZ, and 77 are somewhat smaller than 
those for the SM, these clean signals can be searched for in the near future. For example, 
the Hi could have a sufficient coupling with vector boson pairs to be responsible for the 
approximately 98 GeV excess at LEP [28, 60, 61]. 

Again, we find it very interesting that a heavier Higgs state could dominantly decay to 
a pair of lighter Higgs bosons. Note that Hi is non-SM-like and light, so that there are no 
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Figure 19. Decay branching fractions for Hi 2^2,-, ^1,2 and to the SM particles (and H^, 

in the case of i^2-126. The yehow hnes indicate the corresponding values with the SM couplings. 

Higgs pair channels for it to decay to. We see, from Fig. 20, 



H2 — > HiHi, 






(6.7) 


H3 — > HiHi, H1H2, ZAi, 


AiAi, 


H2H2, 


(6.8) 




H2, 




(6.9) 


Ai^ZHi, ZH2, 






(6.10) 


A2^ZHi, AiHi, A1H2, 


ZH2, 


W^H^, ZH3, A1H3, 


(6.11) 
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(d) (e) 

Figure 20. Decay branching fractions for i^i,2,3, ^1,2 and H"^ to Higgs bosons in the case of 
H2-I26. 



again roughly according to the sizes of the branching fractions at the low values of the mass. 
The collider signatures would be multiple heavy quarks, r's, and multiple gauge bosons as 
commented in the last section. The Higgs pair final states from the decay may serve as an 
important window for a new discovery. 

It was previously noted [16, 62] that in the low-m^ region, the direct production of 
the Higgs boson pairs may be quite accessible at the LHC due to the model-independent 
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gauge couplings for H^H~^ and H^AW^. Additional studies include processes such as 
H2H1 [63], low mass with light Ai [64], two low mass Higgs scenarios [65] and 
Higgs boson pair productions [66]. 

7 Summary and Conclusions 

In the framework of the Next to Minimal Supersymmetric Standard Model, we study the 
Higgs sector in light of the discovery of the SM-like Higgs boson at the LHC. We pay 
particular attention to the light Higgs states in the case when the parameter ruA ^ 2m z - 
Our results, coming from a broad parameter scan after implementing the current collider 
constraints from Higgs physics, lead to the following findings: 

• The Higgs bosons in the NMSSM, namely three CP-even states, two CP-odd states, 
and two charged Higgs states, could all be rather light, near or below the electroweak 
scale (i;), although the singlet-like state can be heavier. The SM-like Higgs boson 
could be either the lightest CP-even scalar as in Fig. 3(a), or the second lightest 
CP-even scalar as in Fig. 3(c), but is unlikely to be the heaviest scalar as in Fig. 3(e). 

• If we relax the perturbativity requirement by allowing the NMSSM parameters A and 
Hi to be larger (see Tables 2 and 3), the allowed region for the mass parameters would 
be enlarged significantly (e.g., black versus green, red and magenta points in Fig. 4, 
Fig. 5 and Fig. 11, etc.). 

• The SM-like Higgs signal at the LHC may be appreciably modified, as seen in Figs. 15(a) 
and 18(a) for production, and Figs. 16(a) and 19(b) for decay. 

• Consequently, the 77 rate can be enhanced (Figs. 6, 8 and 12). The naive correlations 
of 77/yy and VV/bb ratios can be violated (Figs. 7 and 13). Furthermore, if the 
SM-like Higgs can decay to a pair of lighter Higgs bosons, the anti-correlation in the 
VV/bb ratio can be further broken (magenta regions of Figs. 7(b) and 13(b)). 

• New Higgs bosons beyond the SM may be readily produced at the LHC. The produc- 
tion cross sections via gg fusion and VBF could be of the same orders of magnitude as 
those of the SM productions (Figs. 15 and 18). Their decay branching fractions to the 
SM particles could be even larger than those of the SM (Figs. 16 and 19), depending 
on tan/3 and the size of their SM-like Higgs fractions (Figs. 9 and 14). 

• The unique channels for the heavy Higgs signal are the decays to a pair of light Higgs 
bosons (Figs. 17 and 20). The striking signals will be multiple heavy quarks (t, b) 
and tau-leptons in the final states. 

In conclusion, the Higgs sector in the NMSSM provides a well motivated theoretical 
framework consistent with the Higgs boson discovery and the searches at the LHC. The low- 
ruA parameter region yields multiple light Higgs bosons that lead to rich phenomenology at 
colliders. Although strongly constrained by the current searches, it is also highly predictive. 
Dedicated studies for this very interesting sector, in particular for a multiple Higgs boson 
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final state at the LHC, will allow a search for this scenario to be completed in the near 
future. 
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